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We show theoretically that the dynamics of cold atoms in the lowest energy band of a stationary 
optical lattice can be transformed and controlled by a second, weaker, periodic potential moving 
at a constant speed along the axis of the stationary lattice. The atom trajectories exhibit complex 
behavior, which depends sensitively on the amplitude and speed of the propagating lattice. When 
the speed and amplitude of the moving potential are low, the atoms are dragged through the 
static lattice and perform drifting orbits with frequencies an order of magnitude higher than that 
corresponding to the moving potential. Increasing either the speed or amplitude of the moving 
lattice induces Bloch-like oscillations within the energy band of the static lattice, which exhibit 
complex resonances at critical values of the system parameters. In some cases, a very small change 
in these parameters can reverse the atom's direction of motion. In order to understand these 
dynamics we present an analytical model, which describes the key features of the atom transport 
and also accurately predicts the positions of the resonant features in the atom's phase space. The 
abrupt controllable transitions between dynamical regimes, and the associated set of resonances, 
provide a mechanism for transporting atoms between precise locations in a lattice: as required 
for using cold atoms to simulate condensed matter or as a stepping stone to quantum information 
processing. The system also provides a direct quantum simulator of acoustic waves propagating 
through semiconductor nanostructures in sound analogs of the optical laser (SASER). 

PACS numbers: 37.10.Jk, 05.45.-a, 63.20.kd 



I. INTRODUCTION 

Cold atoms in optical lattices created by two or more 
counter-propagating laser beams provide a clean, flexible, 
physical environment for studying quantum particles in 
periodic potentials. Key advantages of such systems in- 
clude the tunability of the lattice, low dissipation, and 
the ability to directly control and measure the atoms' in- 
ternal state. When the atoms are cold enough for the de 
Broglie wavelength, corresponding to the atoms' center- 
of-mass motion, to span several wells of the lattice po- 
tential, the associated energy eigenvalues form bands of 
allowed energies separated by band gaps. The dynamics 
of atoms in energy bands is of particular interest due to 
the similarities with the behavior of electrons in solids 
[IHB]. For example, atoms confined to a single energy 
band and subject to a constant force have been shown to 
Bragg reflect and perform many Bloch oscillations [THU] : 
a phenomenon well known, but harder to realise, in solid 
state physics [HI [H] . 

Recently, there has been increased experimental inter- 
est in nonlinear phenomena relating to the transport of 
atoms in optical lattices. In particular, it has been shown 
that moving optical lattices can transport cold atoms 
along distances as far as 20 cm [T5] . In addition, unusual 
band structure properties have been revealed by experi- 
ments with a moving optical lattice |16| , critical velocities 
in such systems have been identified [17 , the stability of 
superfluid currents in systems of ultracold bosons have 
been studied [18] , and the phenomenon of stochastic res- 
onances has been predicted to play an important role in 



transport through dissipative optical lattices [TO]. For 
single-atom transport, both lensing ^20) and the dynam- 
ics and control [2THS5] of BECs have been observed and 
studied. 

In this paper, we propose a method to control precisely 
the dynamics of atoms moving through a stationary op- 
tical lattice by introducing an additional moving peri- 
odic optical potential. We demonstrate that the com- 
bined effect of the two lattices induces multiple sharp 
resonances, which can be exploited to tailor the atomic 
motion. Also, this system has a direct counterpart in 
semiconductor physics. Specifically, it models the elec- 
tron dynamics induced by an acoustic wave propagating 
through semiconductor nanostructures |27| : a topic of 
growing interest due to the recent development of co- 
herent sources of phonons (SASERs) [251130) . which are 
analogs of the optical laser. 

The paper has the following structure. In Section |llj 
we introduce the semiclassical equations of motion for an 
atom in a stationary optical lattice driven by a propa- 
gating optical potential. Analysis of the equations re- 
veals the distinct dynamical regimes of this system and 
the abrupt transitions between them. In Section [III[ we 
study how the speed of the propagating potential affects 
the average velocity of atoms in the stationary lattice. 
Section [TV| explores how the atom's initial position affects 
transitions between the various dynamical regimes. In 
Section [Vj we compare our semiclassical results with full 
quantum calculations. Finally, in Section |VH we draw 
conclusions and propose experiments to study the trans- 
port phenomena that we have identified. 
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FIG. 1: (a) solid curve shows the potential energy, Vol, of a 
■^^Na atom in the stationary optical lattice. Shaded regions 
show energy ranges of first and second energy bands, (b) 
E(px) curve calculated for the first energy band, (c) Solid 
curve: total potential, Vol+Vm, of the atom ait — 0. Dashed 
curve: total potential after half a temporal period of Vm (i.e. 
at t = -K /ujm)- 



II. SEMICLASSICAL ATOM TRAJECTORIES 

We consider a cloud of cold, non-interacting, sodium 
atoms [m I5TI - I51] initially in the lowest energy band of a 
ID stationary optical lattice (SOL), taking the potential 
energy profile of each atom to be Vq l {x) — Vq sin^ {ttx /d) 
[Fig. [TJa)]. The lattice period d = 294.5 nm and its 
depth is Vq = 563 peV = 5.5^;^, where Er = 103 peV 
is the recoil energy. The energy versus crystal momen- 
tum dispersion relation for the lowest energy band is, to 
good approximation, E{px) = Eo + A[l — cos{pxd/h)] /2 
[Fig. [ijb)], where Eq « 200 peV and A = 24.4 peV 
is the band width. To this system we apply an ad- 
ditional moving optical potential (MOP), which propa- 
gates along the x axis at velocity, vm, creating a position 
and time (t) dependent potential energy field, VMix, t) = 
Um [1 — sin {kMX — uiMt)] /2 [Fig. [T|c)], where the wave 
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FIG. 2: Solid curve shows {vx)t versus the moving optical 
potential amplitude, Um- Vertical (red) dashed line marks the 
onset of atom dragging (labelled Um)- Dot-dashed vertical 
line (blue) marks the onset of Bloch oscillations (labelled Um)- 
Arrows highlight resonant peaks in {vx)t, explained in the 
text. For guidance, the upper dotted line shows {vx}t ~ vm 
and the lower dotted line shows {vx)t = 0. 



number k^i = ^tt/Xm, Am is the wavelength, and lom = 
kiiiVM is the angular frequency. Such a field can be gener- 
ated by counter-propagating laser beams whose frequen- 
cies are slightly detuned by Su- The velocity of the MOP 
is then VM = Suj/2kM fHHHIIl. 

The semiclassical Hamiltonian for an atom in the low- 
est energy band of the stationary optical lattice, and sub- 
ject to the moving potential, is H{x,px,t) = E{px) + 
VM{x,t). The corresponding Hamilton's equations of 
motion are 



dpx 
dt 



U. 



M 



kM COS [knix + xo) - i^Ait] 



dx 

Ik 



Ad^ 
'2h' 



Pxd 

h 



(1) 

(2) 



We solve Eqs. ([I]) and ^ numerically to calculate the 
trajectories starting from rest with xq = x{t = 0) = 
and Vx{t = 0) = 0. To characterise the transport of an 
atom through the SOL, we calculate the time averaged 
velocity, {vx)t, over a period of 0.25 s. 

The solid curve in Fig. [2] shows {vx)t calculated as a 
function of Um for \m = 20d and vm = 2.5 mm . The 
curve reveals that for very low Um, {vx)t increases expo- 
nentially from with increasing Um, until Um reaches a 
critical value, Ufj « 2.5 peV = 2.43 x 10~^Er (red ver- 
tical dashed line), at which point {vx)t = vm (marked by 
upper horizontal dotted line). As Um increases beyond 
^M' {^x)t remains pinned at vm until U]\i reaches a sec- 
ond critical value, C/^ ~ 16 peV = 0.156i?fl:, (blue ver- 
tical dot-dashed line) where the average velocity begins 
to decrease abruptly with increasing Um- Thereafter, 
increasing Um gives rise to a series of resonant peaks 
(arrowed in Fig. [2]). 

To understand the form of the {vx)t versus Um curve. 
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in Fig. [3] we show atom trajectories for three values of 
Um corresponding to three distinct dynamical regimes. 

Figure [Sf^a) shows x versus t for a trajectory within 
the low field regime when Um — 1 peV — 9.74 x 
(< Uff, which is marked by red dashed line in Fig. [2|. 
In this case, the atom orbit comprises a uniform drift, 
with mean velocity w 75 /.tm s~^, and superimposed pe- 
riodic oscillations. The electron velocity > when 
VMix,t) < Um/'^ [white regions in Fig. |3][a)] and t^a; < 
when VM{x,t) > UmI'^ [gray regions in Fig. |3][a)]. We 
explain the form of this trajectory by considering atom 
motion in the rest frame of the MOP, in which the atom's 
position is x' (t) = x{t) — VMt. In this frame, the semiclas- 
sical Hamiltonian is H' {x',px) = E'{px) + Vm{x'), where 
E'{Px) = E{px) — vmPx is the effective dispersion rela- 
tion and the potential, Vm{x') = Um [1 — sin {kMx')] /2, 
which does not explicitly depend on time. The trans- 
formed semiclassical equations of motion are 



dpx 

at 

VM = - 



Um 



kM cos[kMix' + xq)] (3) 



dt 



Ad ^ 



Pxd 

h 



VM- 



(4) 



Since H' is time- independent, it is a constant of motion 
but does not equal the total system energy, H. For initial 
conditions x{t = 0) = xq and Px{t — 0) ^ Pq we obtain 



H'{x',px) 



Um 



[1 - sin {kMXo}] 



(5) 



Therefore, if the initial position of the atom is = 0, 
the constant of motion is H' — Um/"^, which means 
that the kinetic energy in the moving frame E'(px) = 
H'{x' ,px) — Vm{x') = Um/'2' — Vm{x') can only take val- 
ues between ±[/jif/2, marked by the horizontal dashed 
lines in Fig. [3j^b) for Um = 1 peV. This figure reveals 
that if po ~ 0, the atom can only access a very limited, 
almost linear, region of E' [px) near Px = Q (between the 
horizontal dashed hues), where v'^ = dE'{px)/dpx ~ —1^, 
with v being a positive constant. Consequently, in the 
rest frame, the mean velocity of the atom, {vx)t ~ ^m — v^ 
is smaller than vm- This can be seen from Fig. |3][a), 
where the mean (drift) velocity of the atom orbit, is less 
than the slope of the grey and white stripes, which equals 
Vm- As the atom traverses successive minima and max- 
ima [white and gray stripes in Fig. [3](a)] of the mov- 
ing lattice, the corresponding force changes sign, thereby 
causing the atom to oscillate between the extremal px 
values marked by the circles in Fig. pib) and, in real 
space, oscillate around the mean drift [Fig. |3];a)]. As 
Um increases [i.e. the horizontal dashed lines in Fig. 
[3f^b) move further apart] , the atom can access more of the 
E'{px) dispersion curve. As a result, v becomes smaller 
and {vx)t increases towards vm, as shown by the section 
of the black curve in Fig. [2] to the left of the vertical 
red dashed line. Atom motion in this 'Linear Dispersion' 
(LD) regime is described in detail in Appendix [A| 




FIG. 3: x{t) [(a), (c) and (e)] and E'{p^) [(b), (d) and (f)] 
curves calculated for Um = 1 peV = 9.74 x 10~^Er [(a) and 
(b)], 10 peV = 9.74 x W'^Er [(c) and (d)], and 150 peV = 
1A6Er [(e) and (f)]. In (a), (c) and (e), VM{x,t) is < Um /2 
[> Um/'S], within the white [gray] regions. In (b), (d) and (f), 
horizontal dashed lines show ±Um/'2 and filled circles mark 
when E'{px) = ±Um /'2 and, therefore, the limits of the orbits 
in p^. Arrows, labels, and brackets are discussed in the text 
and Appendices [A] and [B] 



The LD regime persists with increasing Um until the 
value of —Um/'^ [lower dashed line in Fig. [sj^b)] falls 
below the value of the local minimum in E'^{px) [ar- 
rowed in Fig. |3][b)]. Using the fact that in this regime 
\Px \ ^ nh/d, we find that the Px value corresponding to 
the local minimum of E'{px) is Px ~ 2T?vm I A.tP' . Sub- 
stituting this value into the expression for E'{px), we 
obtain the following estimate for the critical amplitude 
of the moving lattice, Ufj, above which the atom can 
reach and traverse the local minimum in E'{px)'- 



Um = X 

1 + smkMXo 



(6) 



Ad2 



-A I 



2hv 



M 



Ad 



For the parameters considered, Eq. ^ gives Ufj ~ 2.5 
peV — 2.43 X IQ^^En, which compares well with the 
transition to the region where {vx)t = vm in Fig. |2] (i.e. 
to the right of the vertical red dashed line). 

The solid curve in Fig. [3f^c) shows the trajectory of the 
atom when Um = 10 peV = 9.74 x IQ-'^Er > Ufj. We 
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find regular, almost periodic, oscillations superimposed 
on a uniform drift, indicated by the dashed line in Fig. 
[Sjjc). The slope of the line is wm, indicating that the 
atom is dragged through the SOL by the MOP. This is 
confirmed by the fact that the trajectory is confined to a 
single minimum in the MOP potential where VM{x,t) < 
[gray region in Fig. 3]Jc)], which traps the atom and 
transports it through the SOL. Fig. [sjjd) shows that the 
atom is able to oscillate in a parabolic region of E'{px) 
between the two Px values (filled circles) where E'{px) = 
Um/'^- Since the atom remains in the almost parabolic 
region of E'{px), x'(t) is an almost harmonic function of 
t. Therefore we can approximate the trajectory in Fig. 
[Sjjc) by the expression x{t) « vmI + Xm [1 — cos (ojfli)] /4, 
where uji is the frequency for motion to and fro across 
the potential well In this 'Wave Dragging' (WD) 

regime, the atom is trapped in a single well of the MOP, 
but as the well moves, it drags the atom through the SOL 
with {vx)t = vm (upper horizontal dotted line in Fig. [2]). 

Increasing C/m in the WD regime initially has no qual- 
itative effect on the trajectories, which continue to be 
dragged and are of the form x{t) = vm^ + /(*)■ How- 
ever, as Um increases, the atom can access increasingly 
non-parabohc parts of E'{px) and thus f{t) becomes less 
harmonic. Eventually, though, Um becomes large enough 
for the atom to traverse the first local maximum of E' {pj.) 
[arrowed in Fig. [sj^d)]. Thereafter, the atom can reach 
the edge of the first Brillouin zone and its trajectory in p^ 
changes abruptly from closed to open. In this regime, the 
atom can traverse several minizones, allowing it to Bragg 
reflect and perform Bloch oscillations. The first local 
maximum of ii^' (pa; ) occurs when fm / d—ih^VM / ^d? 
[27] . Substituting this value into E'ip^) we find the fol- 
lowing approximate expression for the wave amplitude, 
C/^, corresponding to the transition between the WD 
and 'Bloch Oscillation' (BO) regimes, i.e. for which the 
value of C/a//2, shown by the upper dashed line in Fig. 
[sjjd) reaches the local maximum in E'{px) P7] : 



M 



1 — sin /ca/Xq 



A 



hrrvM 



Ad2 



(7) 



The vertical blue dot-dashed line in Fig. [2] shows the 
value of U'^ obtained from Eq. ([t]) for xq = 0, which 
coincides exactly with the abrupt suppression of {vx)t- 

Figure [sje) shows the atom's trajectory when Um = 
150 peV = lAGEfj > U\j. In this regime, the atom 
undergoes fast oscillations [bracketed in Fig|3][e)], inter- 
rupted by jumps in the orbit (arrowed). Initially, for 
t = and x{t = Q) = px{t = Q) = 0, the driving force on 
the atom is maximal, and therefore Px increases rapidly. 
As a result, the atom traverses several local maxima and 
minima in E'{px) [between circles in Fig. Islf)], Bragg 
reflecting and performing Bloch oscillations [within left- 
most bracketed region of trajectory in Fig. [sjje)] as it 
does so. As E'{px) decreases, Vm increases to keep H' 
constant. Eventually, Vm attains its maximum value of 
UmI^ corresponding to the minimum attainable value of 



E'{px) ~ —Um^^ [lower dashed line in Fig. Sj^f)] and 
maximum attainable value of Px [right-hand flUed circle 
in Fig. [Sj^f)]. Thereafter, px starts to decrease, trigger- 
ing another burst of Bloch oscillations, until px reaches 
its minimum possible value [left-hand filled circle in Fig. 
[Sj^f)]. Subsequently, Px increases again and the cycle re- 
peats. The details of the form of this orbit are given in 
Appendix |B] 

The analysis in the appendix shows that the atom's 
average velocity will exhibit a resonant peak whenever 
Um just exceeds a local maximum, or falls below a local 
minimum, in E'{px), thus maximising the time the atom 
spends in a region of E' (jpx) where v'^ = dE'{px)/ dpx > 0, 
and so also maximising {vx)t — Wx)t +vm ■ The minima 
and maxima in the dispersion curve occur, respectively, 
at Px values given by 



and 



Px 



<'2r 



2TTh 2h^VM 

'^ + ^d^' 



7r/i 2h^VM 



(8) 



(9) 



where rmin and rmax are integers {— 1,2,3...) labelling 
the minima and maxima, respectively, in E'{px). The 
corresponding critical values of Um, obtained from set- 
ting E'{px) = E{px) - vmPx = ±Um/'2., are 



U 



and 



M 



^ M 



-2E{pl- 



2e{p: 



2vmpI" 



2vmpI" 



(10) 



(11) 



The values of f/]J}°^ and for r^ax 



and Tn 



and 2 are shown in Fig. [2] by the labelled arrows and 
correspond very well with the positions of the resonant 
peaks in {vx)t- 



III. VARIATION OF MOVING LATTICE SPEED 

The parameters of a moving optical potential can be 
varied over a wide range [15l |33]. For example, it has 
been shown that vm can be increased to ~ 50 mm s~^ 
[H] and Um to ~ 500 peV [33]. 

The color map in Fig. |4| shows the variation of the 
average velocity of the atom, {vx)t, with both vm and 
Um ■ The figure reveals a complex patchwork of distinct 
dynamical regimes. At the boundaries between these 
regimes, the average velocity changes abruptly and may 
even change sign, indicating that the atom's response to 
the moving wave depends critically on both the ampli- 
tude and velocity of that wave. 
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FIG. 4: Color map (scale right) showing the variation of {vx)t 
with vm and Um- Dotted black curve shows Ufj versus vm, 
and marks the onset of dragged atom trajectories [Eq. (|6|]. 
Black solid curve shows U\j versus vm [Eq. ij?])]. Black dot- 
dashed curves show (7]^™'" versus vm when rmin = 1,2 and 
3 (bottom to top) [Eq. ( |10[ )]. Red dot-dashed curves mark 
values of when rmax ~ 1,2 and 3 (bottom to top) 

[Eq. (111]. Labels LD, WD and BO mark regions of 'Linear 
Dispersion', 'Wave Dragging' and 'Bloch Oscillations'. 



The structure of the color map in Fig. |4] can be un- 
derstood by considering the transitions between the LD, 
WD and BO regimes discussed in the previous section. 
The dotted black curve in Fig. |4] shows Uf,j versus vm 
[Eq. ^] , marking the transition between the LD and the 
WD regimes. For vm 3.5 mm s~^ (i.e. to the left of 
the intersect between the solid and dotted black curves), 
increasing Um across the black dotted curve induces a 
clear transition from the LD regime, where {vx)t ~ 
(blue in Fig. [4]), to the WD regime, where {vx)t ~ vai 
(yellow / red in Fig. |4]). 

The black solid curve in Fig. ^ shows U'^ versus vj^i 
[Eq. ([7|] and thus marks the transition from the WD 
regime into the BO regime. The region of the color map 
below the black solid curve and above the black dotted 
curve corresponds to the WD regime where {vx)t — vm- 
At the lower edge of this region, the color map changes 
from blue (low {vx)t) to red (high {vx)t) as Um increases. 

As Vm increases, the range of Um over which dragging 
occurs decreases and vanishes at vm = lS.d/hTi = 3.5 mm 
s~^, when Ufj — U\j (crossing of solid and dotted black 
curves in Fig. At this vm value, the tilt of E'{px) is 
large enough to ensure that the value of Um needed for 
the atom to traverse the first local minimum in E'{px) 
is the same, or larger, as that required to traverse the 
first local maximum. Therefore, when Um ~ U'^j, the 



dynamics change straight from the LD to the BO regime, 
i.e. without traversing the WD region. 

In the BO regime, we expect the atom to have its high- 
est value of {vx)t when Um — Ufj since here the atom 
spends most time in the region of the E' {p^) curve where 
the gradient of E'{px), and thus v'^., is positive. This is 
confirmed by the color map in Fig. |4j which reveals a nar- 
row region of high {vx)t values (red), where {vx)t vm, 
just above the dotted black curve along which Um = Ufj. 

Also included in Fig. |4] are black and red dot-dashed 
curves showing, respectively, C/]^"*" [Eq. ( 10 )] and U^]"''' 
[Eq. (11)] for 

fmin and T^ax — 1,2 and 3. These curves 
show good correspondence with the positions of {vx)t 
peaks in the color map, which appear as blue / green 
or yellow stripes for Um > 40 peV = 0.390£'fl;. 



IV. EFFECT OF INITIAL POSITION 

Equations ([6| and ([7| indicate that the transitions be- 
tween distinct dynamical regimes depend on the initial 
position of the atom x{t = 0) = xq- To illustrate this. 
Fig. [5] shows a color map of {vx)t versus xq and U m when 
Vm = 2.5 mm s~^ and, as before, Xm — 20rf. The color 
map has an intricate structure which, we now explain, 
is associated with the transition between the various dy- 
namical regimes. 

The solid curve in Fig. [2] corresponds to a vertical slice 
through the color map in Fig. [5] when xq — 0, i.e. along 
the middle of the five vertical black lines. Increasing Um 
along this line, induces a transition from the LD into 
the WD regime on crossing the black dotted curve where 
{vx)t — Vm- The WD regime appears as the dark red 
region in the color map. With further increase of Umi we 
enter the BO regime (on crossing the black solid curve), 
where the color map changes abruptly from red to yellow 
and the atom's velocity falls dramatically. Within the 
BO regime there is a series of resonances, which appear 
as yellow areas in Fig. [5] and correspond to those arrowed 
in Fig. [2] As explained in Section |ll] and Appendix |B] 
these resonances occur whenever Uj^i coincides with a 
local extremum in E'{px). 

The variation of Uf^ with xq [Eq. (j6])], shown by the 
dotted curve in Fig. [5j agrees well with the transition 
from the LD to the WD regime seen in the color map. 
As Xq increases from to Am/4 the transition to the 
WD regime shifts slowly to lower Ufj, which attains a 
minimum value of 1.3 peV = 1.27 x IQ^^Er when xq ~ 
Am/4. For Xq > Am/4, Ufj increases slowly until xq = 
Am/2, at which point the dynamics are approximately 
equivalent to those at xq — 0. Conversely, decreasing xq 
below makes Ufj increase rapidly until xq ~ — 3Aa//16 
when the dotted curve in Fig. [5]diverges. As Xq decreases 
further, we enter a regime in which there is no transition 
to the WD regime. When xq « — 5Am/16, the dotted 
curve re-emerges and the value of Uf^ decreases rapidly 
with decreasing xq until xq = —Am/2, which is exactly 
equivalent to xq = Xm/2, meaning that the color maps 
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FIG. 5; Color map (scale right) showing the variation of the 
average atom velocity with xq and Um- Black dotted curve 
shows the value of Ufi{xo) representing the onset of dragged 
atom trajectoriesJEq. ([6|]. Black solid curve marks the values 
oiUli{xo) [Eq. |7|]. Black dot-dashed curves show (7]J^™'" (a;o) 
when Tmin ~ 1,2 and 3 [Eq. (|10[)]. Red dot-dashed curves 



are the values of (7]^}"°°" (a;o) where Tmax = 1,2 and 3 [Eq. 



(11 1]. For guidance, vertical black lines show the positions of 
the center of five wells in the SOL, located at —2d, —d, 0, 
d and 2d. Labels LD, WD and BO mark regions of 'Linear 
Dispersion', 'Wave Dragging' and 'Bloch Oscillations'. 



on the left- and right-hand edges of Fig. [5] are identical. 

Equation Q shows that as vm increases from 0, 
Uf.jixo — 0) also increases and, therefore, the range 
of Um values corresponding to the WD regime is re- 
duced. Figure [s] and Eq. ^ show, in addition, that as 
xq — — Am/4, Uf,j{xQ) — oo, meaning that for any given 
Vm and Um the atom will never enter the WD regime. 
We understand this in the effective dispersion curve pic- 
ture by noting, from Eq. ([s]), that when xq = — Am/4, 
H' = Um and so E'{px) can only take values between 
and Um- Therefore, the atom cannot traverse the local 
minimum in E' {p^) when « [arrowed in Fig. [sj^b)] 
and thus cannot enter the WD regime. 

The variation of U\[ with xq [Eq. ([t])], shown by the 
black solid curve in Fig. [s] is in excellent agreement with 
the transition from the WD regime (red region in Fig. [5]) 
to the BO regime (above both the black dotted and solid 
curves in Fig. [s]). As xq increases from 0, U\i increases 
until xq « 3Am/16 where, for the range of Um shown, 
there is no transition to the BO regime above the red 
region. Equation ([t]) shows that as xq — > Am/4, U%[ — >■ 
oo, implying that trajectories with xq = Am/4 are never 
able to Bloch oscillate for any vm or Um- Physically this 
is because when xq = Am/4, according to Eq. E'{p.j.) 
only takes values between —Um and 0, so the atom is 



unable to traverse the local maximum in E' [px) [arrowed 
in Fig. |3][d)] and thus cannot Bloch oscillate. Figure [s] 
reveals that the Bloch oscillation regime reappears in the 
parameter space when xq ~ 5Aa//16. Thereafter, U\j 
(right-hand segment of the black solid curve) decreases 
with increasing xq until xq — Am/2, which is equivalent 
to Xq = —\m/'2' (color maps are identical on left- and 
right-hand edges of Fig. [5| . 

As Xq decreases from 0, the transition from the WD to 
the BO regime occurs at slowly decreasing Um values un- 
til Uf,j{xo) = UIi{xq) (where the black dotted and solid 
curves cross in Fig. [5]). As xq decreases further, the two 
curves cross again when xq ~ — 3Am/8. For xq values be- 
tween these two crossing points, there is no wave dragging 
regime and, for given xq, increasing Um induces a transi- 
tion straight from the LD to the BO regime upon crossing 
the black solid curve. As xq approaches — AAf/4 then, for 
any given Um, the atom cannot traverse the local mini- 
mum in E' {px) [arrowed in Fig. [sj^b)] and thus also does 
not reach the first local maximum in E'{px) [arrowed in 
Fig. [sjjd)]. Therefore, the atom will only change direc- 
tion when it reaches the peak in E'{px) occurring when 
Px —hn/d, resulting in the dramatic change from neg- 
ative to positive velocity [transition from blue to yellow 
regions m Fig. [5] when Um « 45 peV = 0A38Eb- 

As discussed previously, the strong resonant features 
in {vx)t versus Um and xq in Fig. [s] occur when the 



atom accesses new peaks in E'{px) [see Eqs. (10) and 
(llll)] . The variation of UIT'" and UIT"^" with xo shown 



by the black and red dot-dashed curves respectively in 
Fig. [5] for (bottom to top) Tmin and r^ax = 1,2 and 3 
agree well with the yellow regions of the color map where 
{vx)t increases abruptly. Note that near these resonances, 
the color map changes abruptly between yellow and blue 
indicating that {vx)t changes suddenly from positive to 
negative: i.e small changes in xq or Um switch the atom's 
direction of travel. 



V. QUANTUM CALCULATIONS 

In this section, we investigate the full quantum- 
mechanical description of the semiclassical dynami- 
cal regimes revealed in Sections [Tl|lV| The quantum- 
mechanical Hamiltonian for the system is 



92 



2ma dx^ 



+ VoL{x) + VM{x,t), (12) 



where nia is the mass of a single ^'^Na atom. We solved 
the corresponding time-dependent Schrodinger equation 



ih 



dtjj{x) 
dt 



'H{x,t)'ij;{x), 



(13) 



numerically using the Crank-Nicolson method [35_^ . Note 
that we consider the mean-field interaction to be negli- 
gible compared with the external potential. This can be 
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realised in experiment by using an applied magnetic field 
to tune the Feshbach resonance so that the inter-atomic 
scattering length is zero. By only considering a non- 
interacting atom cloud we can present a full and detailed 
description of the fundamental dynamical phenomena ex- 
hibited by the system. If atom interactions are included 
we expect more complex dynamics which would deserve 
further study and analysis. 

First, we found the form of the initial wavepacket by 
evolving Eq. ( 13 1 in imaginary time, using an initial 



Gaussian wavefunction with a full width half maximum 
value, fx = 2c?, to produce a state similar to the ground 
state of an atom cloud held in the SOL by a harmonic 
trap. To study the atomic dynamics, we then integrated 
Eq. ( 13 ) numerically, using this wavepacket as the ini- 
tial state, which, when — 2d, extends across ~ 5 
wells of the SOL. Figure [6] shows the time-evolution of 
the wavepacket in each of the three distinct dynamical 
regimes corresponding to the semiclassical trajectories 
shown in Fig. |3](wAf — 2.5 mm s^^. Am = 20d) (see also 
Supplemental Information 36j ) ■ The central black curve 
in Fig. [6f^a) corresponds to the semiclassical trajectory 
in the LD regime {Um — 1 peV, a;o — 0) shown in Fig. 
[3f^a) and Supplemental Information, Movie 1 Also 
shown are semiclassical trajectories calculated for (from 
bottom to top) Xq = —2d,—d,0,d and 2d, i.e. for orbits 
starting at the centers of the 5 wells that are spanned 
by the initial wavefunction. In this case, most of the 
wavepacket, shown by the gray-scale map in Fig. |6][a), 
does not follow the semiclassical trajectories. Instead, it 
spreads rapidly with only a small fraction of the wave- 
function being dragged through the lattice. 

From Fig. [5] at first sight we might expect that the 
wavepacket lies entirely in the LD regime because, at 
Um = 1 peV = 9.74 x 10~^Er, the vertical fines mark- 
ing the centers of the potential wells spanned by the 
wavepacket all lie in the light green region below the black 
dotted curve, i.e. fully within the LD regime. It is clear 
from Fig. [6]ja), however, that the wavepacket does not 
follow a LD semiclassical trajectory. This is because for 
such a low Um value the spreading of the wavepacket, re- 
sulting from position-momentum uncertainty, dominates 
the dynamics obtained from the semiclassical analysis. 
From the uncertainty principle, fx^Px ~ fi, where Apx 
is the spread of the momentum of the wavepacket, and 
the corresponding zero point energy of the wavepacket 
is E^p = Apl/2ma. When fx = 2d, E^p « 2.5 peV = 
2.43 X 10~^i?fl,. This energy is large enough to allow part 
of the wavepacket to enter the wave dragging regime, as 
observed in Fig. |6|^a) . This zero point energy determines 
a lower limit of J7m below which the semiclassical analysis 
becomes invalid due to being dominated by wavepacket 
expansion. 

Figure [6](b) and Supplemental Information, Movie 2 
|36j , shows the evolution of the wavepacket when Um — 
10 peV = 9.74 X 10~^ii^/j, corresponding to the semiclassi- 
cal trajectory in the WD regime shown in Fig. |3][c). The 
striking similarity between the semiclassical orbits and 




FIG. 6: Gray-scale maps [color online: white = 0, blue = high] 
showing the evolution of the atom density for: (a) Um ~ 1 
peV = 9.74 X IQ-^Er and = 2d (LD regime); (b) Um = 
10 peV = 9.74 X IQ-'^Er and ^ = 2d (WD regime); (c) 
Um = 150 peV = IAGEr and ^ = 2d (BO regime); (d) 
Um ~ 25 peV — 0.243Er, fx — 8d which encompasses both 
the Wave Dragging and Bloch Oscillation regimes. See text 
for further details. Solid curves in (a)-(c) show corresponding 
semiclassical trajectories with, from bottom to top, xo = —2d, 
—d, 0, d and 2d. Movies showing the time evolution of the 
wavepackets are shown in the Supplemental Material |36| . 



the atom density evolution reveals that, in this regime, 
the wavepacket is dragged through the SOL by the MOL. 
The reason for this is clear from Fig. [5] which reveals that 
for Um = 10 peV the full spread of the wave function (be- 
tween the left- and right-most dashed vertical lines in Fig. 
[5]) lies within the red wave dragging regime. For this Um 
value, the moving lattice is strong enough to dominate 
any spreading caused by the zero point energy. 

Figure |6jc) and Supplemental Information, Movie 3 
[55] . shows the wavepacket evolution for Um = 150 peV 
= lA6Eii, corresponding to the BO regime [e.g. the 
orbit in Fig. Isfe)]. The five semiclassical trajectories 
shown in Fig. T6[c) broadly encompass the evolution of 
the negative-going part of the wavepacket. One might 
expect that the energy band would break for this large 
value of Um, causing the wavepacket to no longer track 
the semiclassical dynamics. However, we note that the 
relatively long wavelength A7\/ = 20d of the moving lat- 
tice ensures that the energy drop across a single well 
within the SOL is « 45 peV = OASSEr, which is small 
enough to keep the band intact. Note also for high val- 
ues of Um the approximation that the atom will remain 
in the lowest energy miniband becomes invalid, we esti- 
mate this upper limit on Um by assuming that interband 
tunnelling will occur when the amplitude of the moving 
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wave is approximately equal to the energy gap between 
the first and second minibands. The two minibands cal- 
culated for the stationary optical lattice are « 250 peV 
apart [see Fig. [ija)]. Therefore we expect interband tun- 
nelling to occur when Um ^ 250 peV. It appears, from 
our numerical simulations, that for Um ^ 200 peV — 
1.95i?ij, the wavepacket dynamics start to deviate from 
the semiclassical paths due to band breakdown. 

The sensitive dependence of the atom dynamics to the 
initial position, xq, shown in Fig. [5j has an interesting 
consequence. Namely, if the atom cloud is large enough, 
it can simultaneously span more than one dynamical 
regime. For example, in Fig. ^d) and Supplemental 
Information, Movie 4 [36], we show the wavepacket evo- 
lution for Um = 25 peV = 0.243-Efl, with an initial 
atom cloud width fx — 8d and center of mass posi- 
tion Xq — l.5d. The figure shows that the atom cloud 
splits, with the lower part performing Bloch oscillations 
around x = and the rest being dragged through the 
SOL by the MOL. The reason is, as can be seen from 
Fig. [5| that the upper part of the initial atom cloud 
with X > 1.5ci lies in the wave dragging regime whilst 
the lower part with a; < 1.5(i is in the Bloch oscillation 
regime. Our semiclassical analysis suggests that approx- 
imately half the cloud is in the wave dragging regime. 
This estimate agrees well with the numerical evolution 
of the wavefunction, which reveals that after 10 ms, the 
fraction of the cloud in the wave dragging regime is 0.495. 
This result suggests that information about the shape 
of the initial atom cloud could be obtained by measur- 
ing the number of atoms, Njj, that end up in the upper 
(dragged) part of the split cloud, as a function of xq. 
For example, as xq increases so that more of the initial 
wavepacket lies within the red WD regime in Fig. [5] i.e. 
so that the initial wavepacket is shifted across the black 
solid curve, more atoms will end-up in the upper part 
of the split cloud. More quantitatively, we expect that 
Nu(xq) « [°°iTT J^(a^o, i = 0)P c?x where xs satisfies 

U\j{xs) — Um, i-e. marks the position of the solid curve 
in Fig. [5] for given Um- 

VI. CONCLUSIONS 

We have shown that a periodic optical potential mov- 
ing through a stationary optical lattice can control the 
atom orbits in real and momentum space, and induce 
transitions between three distinct dynamical regimes. In 
a semiclassical picture, the atom can be confined to a 
small region of the stationary optical lattice dispersion 
curve, dragged through the stationary optical lattice by 
the moving wave, or perform Bloch-like oscillations. The 
crossover between these orbital types creates a rich pat- 
tern of transport regimes, with multiple resonances that 
cause abrupt changes in the magnitude and sign of the 
atom's average velocity. These resonances occur when- 
ever there is an abrupt change in the number of Brillouin 
zones that the atom can access. Since the form of the 



atom trajectories is predictable, and can be tailored by 
small changes in the MOP parameters, the system pro- 
vides a mechanism for moving atom clouds between well- 
defined points in stationary optical lattices. It may also 
provide a sensitive method to split atom clouds, with 
potential uses in atom interferometry and in mapping 
the initial atom cloud dynamically rather than via opti- 
cal imaging. The dynamics described in this paper are 
similar to those expected for electrons in acoustically- 
driven superlattices [27]. Consequently, the cold-atom 
system that we have studied also provides a quantum 
simulator for acoustically-driven semiconductor devices 
in which the electron orbits are much harder to observe 
directly. 



Appendix A: Detailed description of atom motion in 
the Linear Dispersion regime 

In the LD regime, initially, for t = and x' = 0, 
the force applied to the atom by the moving optical lat- 
tice, is maximal, which makes px increase until E' {px) — 
— Um/'^ [lower filled circle in Fig. [sjjb)]. At this point, the 
atom's velocity in the rest frame, Vx — dE' {px)dpx +vm, 
attains its maximum value for the specified initial con- 
dition. The potential energy Vm{x') — Um (to keep H' 
constant) is also maximal, and so the force is zero in the 
moving frame. However, in the moving frame x' contin- 
ues to decrease, which eventually makes the force nega- 
tive {—dVM{x')/dx' < 0), so accelerating the atom in the 
negative Px direction. 

The atom continues to move in the negative px direc- 
tion until E'{px) = Um/'2 [upper filled circle in Fig. [sf^b)] 
where the atom's velocity in the rest frame is minimal. 
The atom then returns to its initial condition with Po = 
and the cycle repeats with the atom traversing successive 
regions where Vm < Um/'2 [white regions in Fig. |3ja)] 
and Vm > Um/'2 [gray regions in Fig. [sj^a)]. Since the 
region of the dispersion curve accessed by the atom is 
not completely linear, the velocities at the extremities of 
the orbit are slightly different, causing the atom to have 
non-zero average velocity. 



Appendix B: Detailed description of atom motion in 
the Bloch Oscillation regime 

In the BO regime, the atom can access E'{px) across 
several Brillouin zones, which makes {vx)t vary in a com- 
plex way with increasing Um- For Um = 150 peV = 
lAGEji, the lowest possible value of E'{px) [lower dashed 
line in Fig. [sjjf)] occurs just below a local minimum in 
the curve [right-hand arrow in Fig. [sjf)]. The velocity 
of the atom in the moving frame, dE'{px)/dpx, is neg- 
ative at the right-hand filled circle, indicating that the 
atom is moving in the negative x' direction. However, 
since Vm^x') is maximal at this point, the moving lat- 
tice exerts little force on the atom. Consequently, px de- 
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creases only slowly and, as a result, the atom also lingers 
in the region of E' [px) where v'^ (and thus Vx) is positive 
[within the region marked by the right-hand bracket in 
Fig. [Sj^f)]. Overall, as the atom moves away from the 
right-hand filled circle in Fig. [sjjf), it spends more time 
with positive than negative velocity causing the atom to 
jump forwards along the x— axis [arrows labelled 'F' in 
Fig-ie)]. 

Conversely, the maximum attainable value of E'(j)x) 
[upper dashed line in Fig. |3]^f)] occurs just below the 
local maximum in E'(px), marked by the left-hand arrow 
in Fig. [sj^f), where Vm is maximal. Since the force on 
the atom is low in this part of the E' {px) curve, the atom 
remains in the negative velocity region of E'ipx) [left- 
hand bracket in Fig. |3]^f)] for a long time. This makes 



the atom take a large jump backwards along the 2;— axis 
[arrows labelled 'B' in Fig. [sjje)]. Figurejsjje) reveals that 
for Um = 150 peV, the atom jumps backwards further 
than it jumps forwards, so giving the atom an overall 
negative average velocity. 

If Um has a value that just exceeds a local maximum or 
falls just below a local minimum in the E' {px) curve, then 
the time that the atom spends in the region of the E'{px) 
curve, where dE' (px) / dpx = v'^ > 0, is maximised. Con- 
sequently, the magnitude of the forward jump in the or- 
bit, and therefore {vx)t = {v'x)t+VM, is maximised. This 
phenomenon produces the series of peaks observed in the 
{vx)t versus Um curve [see arrowed peaks in Fig. [2], the 
positions of which are discussed in the main text. 
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